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Abstract: We prove that the purely classical gravity dual of Fermi and no n- Fermi 
liquids exist by employing the logarithmic behavior of entanglement entropy to probe 
Fermi surfaces. For isotropic systems, the logarithmic behavior originates only from 
the deep UV region of the minimal area surface. For anisotropic systems, the surfaces' 
configuration becomes largely modified by spatial anisotropy and series of solutions 
exist. By imposing the null energy condition we show that the specific heat behaves as 
C oc T a where a < 1, in both systems. In the end, we also present an effective gravity 
model for anisotropic background. 
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1. Introduction 

Recent years, AdS/CFT correspondence has been widely used to study strongly cou- 
pled fermionic systems |], 0, ffl, ]5L 0, [/J. It was shown that holography can describe 
quantum liquids, especially non-Fermi liquids which is still lack of a proper description 
in condensed-matter theory (CMT). By adding probe fermions in the bulk, people can 
explicitly calculate the boundary retarded Green function Gr(oj, k) and analytically 
show the existence of Fermi surfaces [|, || (7j. However, the number of Fermi sur- 
faces are only of order 0(1) in this procedure and usually the Luttinger relation which 
equates the total charge density to the volumes enclosed by the Fermi surfaces was 
badly violated [|], |[, except for the electron star background 0, 0. What's more, 
the bulk fermions only correspond to the gauge-invariant operators in the boundary 
and thus the dual Fermi surfaces are only of these composite operators. Up to now, 
One cann't directly manipulate the Green function and detect the Fermi surfaces of 
boundary elementary fermionic operators in holography |9|, [HJ. Thus we call these 



Fermi surfaces hidden. Interestingly, in a recent paper, N.Ogawa et al have proposed 
an elegant approach to this conundrum . They define systems with Fermi surfaces by 
requiring that their entanglement entropy has the logarithmic violation of the area law. 
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The number of Fermi surfaces is expected of order 0(N 2 ) and the Luttinger relation 
can be partly rescued in some way [|L(| . 



For d spatial dimensional systems, a strip shape subsystem A is defined by: 

i t 

A = {(xi,x 2 , - - < xx < ~,0 < x 2 ,x 3 , ...,x d < L} (1.1) 

When the size I of the subsystem A is large enough, the entanglement entropy Sa will 
be substantially modified and behaves like: 

Sa = + i)L d - x k d i x log{lk F ) + 0(£°) (1.2) 



where 7 and r\ are numerical constants, and rj is positive, e is the UV cut off, \zf is the 
Fermi momentum or the average of Fermi momentums when many Fermi surfaces exist. 
On the other hand, in holography, the entanglement entropy is given by |l^, [13|, [14j 



Area(g A ) 

Sa = 4GV (L3) 

where Gn is the bulk Newton's constant and ja is the minimal area surface which coin- 
cides with dA at the boundary. With the assumption that the logarithmic behavior of 
the entanglement entropy Sa originates from the IR contribution of the minimal area 
surface, Ogawa et al show that only non-Fermi liquids have classical gravity duals. By 
imposing null energy conditions which leads to additional constraints on the metric, 
they show that the specific heat eventually behaves like C oc T a ,with a < 2/3 (d = 
2) and a < 3/5 (d = 3). Thus, only part of non- Fermi liquids can be constructed in 
this way. The standard Fermi liquids (a = 1) and some other non-Fermi liquids are 
not allowed, possibly having no purely classical gravity duals. Is this true? 
Since working in the IR geometry means it is simply to extract the IR pieces of the 
full minimal area surface. It is possible that the other omitted regions is the dominant 
part of the surface area and leads to the logarithmic behavior of the entanglement 



entropy [11]. Motivated by this idea, we will explicitly show that it indeed works suffi- 
ciently for general isotropic systems with d spatial dimension when the deep UV part 
of the surface is dominant (d > 2) in this paper. The index of specific heat will be 
extended to a < 1 and thus include Fermi and all of non-Fermi liquids. 

On the other hand, since the realistic systems in the boundary are generally 



anisotropic, people have recently constructed anisotropic black branes|l7| in holog- 



raphy to further study the properties of the boundary systems. From this perspective, 
it is also important to investigate the hidden Fermi surface information and confirm 
the existence of the purely classical gravity dual of Fermi and non- Fermi liquids in the 
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anisotropic gravity background. We find that the configuration of the minimal surface 
is controlled by two functions product so that many series of solutions which corre- 
sponds to either IR or UV piece of the minimal area surface dominant case are allowed. 
The Fermi and non-Fermi liquids are also naturally allowed in this case. Finally, we 
also present a classical gravity dual for anisotropic systems. 



2. Isotropic systems 

For self consistency, we first briefly demonstrate the calculation process given in Ref||. 
To be general, we consider d + 1 dimensional boundary systems with Fermi surfaces 
which is dual to d + 2 dimensional gravity backgrounds. The metric is taken to be 

ds 2 = ^(-f(r)dt 2 + g(r)dr 2 + dxf) (2.1) 

where R is AdS radius; i = 1, 2, d. Note that this metric preserves translational and 
rotational symmetry and thus is dual to isotropic systems in the boundary. We also 
require this metric is asymptotically AdSd+2, so we have 

/(0) = </(0) = l (2.2) 

We are only interested in the strip subsystem A (1). For other shape subsystems, one 
can refer to Ref|3|, [J. The minimal area surface ja can be specified by the surface 
X\ = x\(r). Without loss of generality, we take the Newton's constant Gn = 1/4 in the 
following. According to eq.(|1.3[), the holographic entanglement entropy is found to be 



-^g^+x'^rf (2.3) 

where r* is the turning point which leads to x[ divergent, e is the UV cut off. Since ja 
is the minimal area surface, the variational principle for x\{r) gives rise to 



*iw = pi/r^ (2-4) 

Thus, we obtain the width I and the entanglement entropy Sa as follows 

(2.5) 





SU = 2«-W~M (2.6) 
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Obviously, the relation between the entanglement entropy Sa and the width £ is con- 
trolled only by the function g(r). Assume that the IR piece of the minimal surface ^a 
is dominant, we can set g(r) as 

g(r) ~ (— ) 2n (r > r F ) 
r F 

~1 (r<r F ) (2.7) 

where r F is a certain scale which is undetermined. Substitute eq. (|2T7j) into eq. (|2.5|) and 
), we find 



„n+l 



(2.8) 



SA-lt-^ + d^L^^— (2.9) 

t I F 

where c n and d n are positive constants and 7^ = 2 d ~ 1 R d /(d — 1). Note that in the very 
IR limit, we have r* 3> r^, or equally £ ^> t>. The eq. (|2.8|) and (|2.9|) lead to 



Jdr-X R d L d ~ l t n-(d-l) 

^ = 7,— + *»-^r - ^ + 0{f) (2.10) 

The first term is divergent and agrees with the area law which is expected from asymp- 
totically AdS background. The second term is the leading finite part of Sa which 
depends on £ as a pow law. An exception occurs when n — d — 1. In this case, the 
behavior of Sa is indeed obtained as 

j d— 1 JDdjd—l p 

Sa = + fc*-i^£r-M-) + O(£ ) (2.11) 

This is the logarithmic behavior we expected for the existence of Fermi surfaces. Com- 
pare eq.( p.llj ) with eq.(2), the scale parameter rp can now be interpreted as the Fermi 
level ~ hp 1 or the average of Fermi levels when many Fermi surfaces exist. 
Up to now, we have only fixed the function g(r) by requiring the logarithmic behavior 
of the finite part of the entanglement entropy while the tt-component of the metric 
Q2.1|) i.e. the function f(r) isn't involved in. To make sure what quantum liquids we 
have in this background, an additional and physically sensible condition is needed. For 
this purpose, let's impose the null energy conditions 

T^N»N U > (2.12) 

where T^ v denotes the energy stress tensor of matter fields; is any null vector. In 
the absence of a specific matter field, T^ u can be calculated from Einstein's equations 

R^u - -g^R = 2vrT^ (2.13) 
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For simplicity, the null vector iV M can be chosen as 

JV = -U J\T = ^, N* = (2.14) 

V 9tt \/9rr 9x^X1 

where 9 is an arbitrary constant. Obviously, equations ( 2.12 )-( 2.14|) apply to both 
isotropic and anisotropic systems. Focus on the IR geometry and assume 

/(r)~r~ 2m , g{r)~r 2n (2.15) 



For isotropic systems (2.1), the above conditions lead to 

m > n and m > (2-16) 
In the end, we found the behavior of the specific heat as [[3] 

C oc T a (2.17) 

m + n + 1 

To show the logarithmic behavior of Sa, we need n = d — 1. Therefore 

d 

a < ^ (2.19) 

Note d > 2. When d = 2, a < 2/3 and d = 3, a < 3/5 which is precisely consistent 
with H]. Clearly, eq.( ^.19[ ) contains only a portion of non- Fermi liquids. The standard 
Landau-Fermi liquids and other non-Fermi liquids are not included in this procedure. 
To generalize above results, one needs to consider the contribution of the deep UV 
region, instead of the deep IR piece, of the minimal area surfaces. We start from a new 
assumption 

9(r) = 1 + zsa ( 2 - 2 °) 



Note that in the deep interior r>rp, eq.( |2.20D reduces to eq.( |2.7] ). Set A = r*/rp ^> 1 
and u = r/r*. Substitute eq.( [2.20|) into eq.(p.5|) and ([2.6|), we get 



. , . „ l + (Au) 2n 
It f = 2 A / duu a - 



u 2d 



~c n A 5 + ... (2.21) 

where c n is a numerical constant, and S is a positive constant related to n in some way. 
The omitted terms is of order 0(A 5_1 ) which seems to be rather bad. Fortunately, 
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log {£/r F ) ~ c„log A + O(A ) which is what we truly need to estimate the finite part 
of the entanglement entropy. In the same way, Sa is obtained as follows 



5A " rt 1 U. 1-^ (2 ' 22) 

2 d-X R d L d-l rl du 



— yi + (A M )2n (2.23) 

T* Je/r* " 

where w 2d term in the integral has been omitted in the second line. Clearly, the m m term 
is dominant when r — > r* which locates in the deep IR region. Since we have assumed 
that the deep UV region contributes to the leading finite part of the entanglement 
entropy Sa, the dominant scale is actually far away from r*. So we expect ( 2.23[ ) is 
sensible. We will show this point explicitly in the following. 

From eq.( 2.23|) and by simple calculation (see Appendix), one find if n = (d— l)/2, we 
will have 

TDd r d— 1 T)d T d—1 
_ it Li . IX Ll . 

S A = ld d _ 1 +d n d-1 logA + ... (2.24) 

which agrees with eq.( 2.11|) perfectly. We emphasize that the dominant scale for this 
behavior to emerge is given by r uv = r F /A s r F , thus locating in the deep UV 
region with the UV cut-off scale e fixed, where 5 is some positive constant to ensure 
this relation. Now let's estimate the leading finite corrections 5(Sa) given by u 2d term. 
When r — > r*, eq. (|2.20|) reduces to eq.( [2.7|) , eq.( p.22| ) reduces to eq. (|2.9|) . Therefore 

pdrd—l , . 

5(S A ) - ^A-(-) (2.25) 

r F 

Thus, for any d (d > 2), 5SA/{SA)fin *C 1, where (SU)/m denotes the finite part of Sa, 
the u 2d term's effect is indeed negligible. 
Combining these results, we finally obtain 

TDdjd— 1 TDd j d—1 p 

Sa = Td^rT- + e «^TT- lo S ~ + ( 2 - 26 ) 

From conditions (|2l2|) - (|2~T8D and now n — (d— l)/2, we find the specific heat behaves 
like 

C oc T a , a < 1 (2.27) 

Thus in this case both the Fermi liquids (a = 1) and all the non-Fermi liquids (a < 1) 
can be constructed. We plot the logarithm of the width log (£/r F ) and the finite part 
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Figure 1: The logarithm of I (log£) as a function of log A. Left plot for d = 2 and right plot 
for d = 3. The blue dotted line shows the numerical result of eq.( 2.21| ), the red line is the 
fitted curve. 
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Figure 2: The finite part of Sa (denoted by (S^/m) as a function of t. Left plot for d = 2 
and right plot for d = 3. The blue dotted line shows the numerical result of eq.( [2.22 ), the red 
line is the fitted curve. 



of the entanglement entropy (SU)/m m Fig. 1 and Fig. 2 1 for d — 2 and d = 3. In 
both cases, we have set rp = 1, R = 1, L = £, A e [50, 100] and rjjv = rp/A. Their 
behaviors are well approximated by 



For d = 2. And 



\og£ 

{S A )fin 

log£ £ 

(5U)/m K 



a 1.52119 logA 

n 103.09£log£-518.124£ 

1.93226 log A 

13838.8£ 2 \ogi- 103972£ 2 



(2.28) 
(2.29) 

(2.30) 
(2.31) 



^^As it was shown in these two plots, there appears a systematic deviation for the numerical points 
from the fitted curve. It behaves so because we drop the small quantity which is of order O(A ) in 



the fitting equations (2.28)-(2.31) 



For d = 3. Thus we can numerically confirm the analytical approach given above is 
indeed rational, as we have expected. 



3. Anisotropic systems 

In the sections above, we discussed isotropic systems with Fermi surfaces in holography 
by requiring the logarithmic violation of the entanglement entropy. When the size of 
the systems is large enough, this behavior is proposed to be general, independent of 
the space shapes of the systems || . In the absence of a proof in the boundary theory, 
we argue that this behavior is also independent of the shape of the momentum space 
i.e. the shape of the Fermi surfaces. In fact, people have found this behavior by simple 
dimensional analysis []15 1. So eq.( |1.2| ) may be also effective for anisotropic systems only 



if we reinterpret the scale parameter kp as the average of the Fermi momentums in 
different directions. 

To put it simply, we are only interested in the two spatial dimensional systems in the 
following. The new metric is taken to be 

ds 2 = ^j(-f(r)dt 2 + g(r)dr 2 + b(r)dx\ + a(r)dx 2 2 ) (3.1) 
We require that this background is also asymptotically AdS, thus 

/(0) = (y(O) = a(0) = b(0) = 1 (3.2) 
And the subsystem A is still a strip one 

A = {(x 1 ,x 2 )\-^<x 1 <^,0<x 2 <L} (3.3) 
By the same procedure presented in section 1, we find 



/ gjr) 

r 2 V a(r)b(r) 2 - b(r)(£) 



*i=W . .„ a:".. (3.4) 



where r* is defined as the turning point of the minimal surface which makes x[ divergent, 
thus we have 

a(r*)6(r*) = 1 (3.5) 
To meet this condition, we first assume 

a{r)b(r) = 1 (3.6) 
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Therefore, we obtain 



r* r 2 ja(r)g(r) 



o 



S A = 2R 2 L 




(3.7) 



(3.8) 



Compare eq.([3.7|) and eq.( |3.8|) with eq.( |2.5| ) and eg. (|2.6| ), we only need to replace g{r) by 
a(r)g(r) when moving from isotropic systems to anisotropic ones. Thus the argument 
presented in section 1 can be directly used to probe the anisotropic Fermi surfaces. 
From eq. fl2.7| ) and eq.( |2.20p , we find that there are several series of solutions to show 
the logarithmic behavior of the entanglement entropy. We list them in the following 



a(r)g(r) 



r F 



a(r)g(r) — 1 + 



IR 



UV 



a(r)g(r) 



'I£)(M-4)r 1 + 
v r 



r 



d>3; UV 



(3.9) 

(3.10) 
(3.11) 



r F 

where IR, UV symbol denotes the cases of which part of ja dominating the area integral. 
Apparently, eq.( [3.9|) and eq. (|3.10|) is from eq.( [2.5|) and eq.( [2.6|) by setting d = 2 while 
eq.( p.ll|) needs further analysis. Substitute eq.( [3.11|) into eq. (|3.8|) and compare with 
eq. ( |2.22|) , the two integral still differ by a functional factor a/(T~ 



u 



2d 



)/(l 



U 



4 ), where 



u = r/r* G [0, 1]. One can prove that this function is monotonic increasing and varys 
from 1 to a/g?/2. For any finite d which is of order 0(A°), this function can be estimated 
by its average which contributes only a constant factor in eq.( |3.8j ). For large d limit, 
this estimation is however broken. But, if the u 2d and u 4 terms can be omitted, we can 
still get eq. (|3.11| ). The point is the dominant contribution to Sa is from the deep UV 
region. 

Furthermore, if we always believe this point, the condition (|3.6|) is no longer needed. 
Thus, we will find 



S A 



2R 2 L 



2R 2 L 




b{r)a?{r)g{r) 



b{r)a{r) 



\/a(r)g(r) 



74 



(3.12) 



Compare eq. (gT2|) with eq.(|2~23|), we again obtain eq.( p,10|) and eq. (|3. ! The big 
difference is we have much more freedom to select proper functions a(r) and b(r) now 
since we only need to meet the weaken condition ( |3.5|) instead of eq.( |3.6|) . For example, 



- 9 - 



we can set a(r)b(r) = 1 only holds in deep IR which is clearly weaken than eq.( [3.6|) but 
makes eq. ([3.5|) satisfied automatically. 

Note that in the case of eg. (|3.9|) and eg. ( |3 . 1 0|) (with condition ( |3.6|) ), the width and the 
entanglement entropy is precisely egual to the one of isotropic systems while in the case 



of eg.( |3.10[ ) (when condition ( |3.6| )) broken) and eg.( |3.11| ) only their logarithm keep the 



same order of A. In form, eg. Q3.llD reduces to eg. ( |3.10| ) when d = 2 and to eg. ( p. lip 



when d = 1 which differ by an unimportant 2 factor. Thus eg. ( 3-11 ) represents the 
general configuration of the minimal area surface ^a- 

To specify the behavior of the specific heat, we impose null conditions on the metric 
(ET1D . From eguations (^T2|)-(^14|), we get 

9' f a' 2 
rg rj 2a z 

(I L^-i — + V— (- >0 (3.14 

K 2f Afg rf Ap } l 2a 2a 2 a l 4/ Ag r )l ~ 1 ' 

We only interested in the deep IR region and assume 

f(r) oc r~ 2m , g(r) oc r 2n , a(r) oc r 2p (3.15) 

which leads to 

m > n + p 2 (3.16) 

(m - p)(m + n + 3) > (3.17) 

On the other hand, according to eguations (|3.9|) - (|3.11|) , the logarithmic behavior of Sa 
reguires in the deep IR 

n + p= (3-d)/2 , d> 1 (3.18) 
For physically sensible matter fields, a > which is eguivalent to 

m + n + l>0 (3.19) 

In the end, we obtain 

2 

d=l, a< — - — - , -oo < p < +oo , m > p 2 - p + 1 (3.20) 

(p -iy + 2 



d = 2, a < 



(p - l) 2 + 1 
a < 4/3 , 1 - y/2/2 <p<l 



p > 1 + V2/2 or p < 1 - V2/2 , m > p 2 - p + 1/2 

2/2 , m > p (3.21) 
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2 2 
a — 3, a < — , p>2orp<0,m>p — p 

(p - iy 

a < 2 , < p < 2 , m > p (3.22) 



d = 4 ' a - (p _ 1)2 _ i ' p - 1 + or p < 1 - ^/372 , m > p 2 - p - 1/2 



a 



< 4 , 1 - y/3/2 <p<l + , m>p (3.23) 



d > 5, a < — , p > p 1+ or p < pj_ , m > p 2 - p + (3 - d)/2 

(p - iy + (3 — d) 

a < +oo , p 2 - < p < P2+ ,m > p + (d - 5)/2 (3.24) 

where p\± = 1 ± y/d — 3 and p 2 ± = 1 ± \/(d— l)/2. From equations ( p.20| )-( ^24] ), one 
can fix a by selecting proper p and m in the permitted interval. Obviously, when d = 1, 
we have a < 1, thus Fermi and no n- Fermi liquids can be constructed but when d > 2, 
even a > 1 is allowed in (|3.21|) -( |3.24|) which represents the thermal systems which we 



have no interests in. Note that every case above suggests that there exists gravity duals 
of Fermi and non- Fermi liquids in the purely classical limit for anisotropic systems. 
On the other hand, above solutions apply to a(r)b(r) = 1 (at least in deep IR) and 
a[r) varies with the scale parameter rp so that in the bulk interior it behaves like 
air) ~ r 2p which is an unnecessary condition. Especially, when the anisotropy is only 
perturbations on the isotropic background, the behavior of a(r) can be set as 

a(r) = 1 + S(a(r)) = 1 + a r 2p (3.25) 

5(a(r))/a(r) < 1 , for all r (3.26) 

where do, p are positive constants. Under this assumption, a"(r)/a(r) ~ a'(r)/a(r) ~ 
5(a(r))/a(r) which can be neglected in the null energy conditions. Thus we obtain 
eq. fl2.16|) again. On the other hand, equations ( |3.9| )-( [3~Tl~1) lead to 



n=(3-d)/2 + 0(6(a)/a) (3.27) 

Hence 

2 

a < - — - + 0(5(a)/a) , d — 1,2,3 , m > n 
a < 4 + 0{5{a)/a) , d = 4 , m > 

a < +oo , d > 5 , m > -{n + 1) (3.28) 

Clearly, for d > 2 the standard Landau-Fermi liquids has been included while for d = 1 
only non-Fermi liquids is allowed. 
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4. Gravity Model 



Now we would like to give an effective gravity model for systems ( |2.1| ) and fl3.ip which 
have the logarithmic behavior of the entanglement entropy. The system ( |2.1| ) is isotropic 
and has been completely constructed by standard Einstein-Maxwell-dilaton theory B 



lq, p77| . So we only focus on the anisotropic system (|3TT| ) with a{r)b{r) = 1. We start 



from the following action 

S = ^J dd+2x ^^ ~ 2 ^<W ~ V i ( t>) ~ y(ZiO)^ 2 + Z 2 { ( j ) )H 2 } (4.1) 

where G denotes the metric determinant, 3ft is Ricci scalar, V((f)) and Z 12 (0) can be 
any function. F,H are two Maxwell fields' strength defined by F = dA, H = dB, 
where A, B are gauge potential. 
The equations of motion(EOM) are 

G> = (4.2) 

V M (Z 1 (0)F^) = (4.3) 

V,(Z 2 (<P)Hn = (4.4) 

where is covariant derivative operator, A is Laplacian operator. G pv is Einstein 
tensor. T pv is energy-momentum tensor which is given by 

Tfiu = 2<9 M <R0 - g, u m) 2 + \v{<t>)) + K 2 Z 1 (4>)(g x »F flX F vp - \g^F 2 ) 

+K 2 Z 2 ((P)(g Xf) H, x H up - \g^H 2 ) (4.6) 



By definition, G^ v of the metric (|3TT| ) are found to be 

r f'(9'^^ M fy a' 2 a> f g> 1., f , . 

°tt = + - + 7T o 2 + _ 77 ~a + 7~ fi r ; 

rg g r 2g 2a a z a 4/ Ag r 2a 

r ^L-°L^ lL<L - I- - \\ + A _ Inn 

r/ 4a 2aHf Ag r } + r z 2a ^ 

1 rr a" a'/' a' 1,. . f' 2 fg' 2f f" 2g' 6,-, 1 0/ , 

G X1X1 = — { — + - 77 "J + "772 - 7f--7- + J y + + ->- ^ r ; 

2ga 2a a 4/ Ag r 2j z 2jg rj j rg r z 2a z 

Gx2X2 ~ 2g i[ 2a a z + a { Af 4<? r )J { 2f z + 2fg + rf f rg r zH + 2 W; 

(4.7) 
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where Q(r) is 

. . a' 2 a' a" a' f a'q' 

tyr) = — + -- ^- 7^-7^; ( 4i 
cr ra 2a Aaj Aag 

To realize the anisotropic metric (ft.l|) , one needs T^ u = T^ir) and T XlXl ^ T, 
T X1X2 — O[0. We search for the following solutions 



X2X2 5 



(j) = <j){r) , A = A(r)dt , B = B(r)dx 1 (4.9) 
EOM ( f4.3|) -( f4~4|) lead to the only nonzero component of gauge fields strength is 



, , H rxi = (4.10) 

V^GZ^) V=GZ 2 (cf>) 

where Qi,Q 2 are constants which are proportional to the corresponding conserved 
charge carrying by the blackholes. For convenience, we set Qi = Q2 — 
Equations flC3)-( |Eiq) lead to 



i'2 



a f a' a' 2 1 

JLQ( r ) _ (J + A_ + " + ^e(r)) (4.11) 

2a V ; V 2/r 2#r 4a 2 2 v ;; v ; 

r 2 f' 2 fV 2f f" a' fi 

= f ( T« + 7^ + 7^-77---^ + ( r ))> ( 412 ) 
2# 4/ 2 4fa /r 2/ gr r 2 

1 f" f' 2 fa' 2f 1 

= 7-^(T - 777 " 7T - T") + — (O; (4-13) 
2gr z j 2j z 2jg jr ar z 

Z-\ r ) = ^O(r) - 4«(r); (4.14) 

0(r) = 2^ + a ( 47-4^-^ (415) 



where 



Note the dilaton equation ( |4.5| ) remained to be a constraint. Generally, it is not satisfied 
automatically. Substituting ( |4.11| )-( fi.l4| ) into it will result to a constraint condition on 
the metric. Hence, we can't have a complete set of solution to all EOM all the time. 
Instead, only a subset solution is allowed. This is quiet different from the case of 
isotropic systems^. To make it clear, let's consider the deep bulk interior: the metric 
components will behave like eq.( |3.15|) , combined with eq. Q4.5D and (|4.11|) - ([4.14|) , we 
obtain 

(2m + 4)p + 3m + n + 9 = (4.16) 

p(2p + m + 3n + 5) - (n + l)(m + n + 3) = (4.17) 

which presents a specific condition on the IR solutions of the ansatz (|3.15|) . Only if 
the index m, n, p locates in the restricted region, our effective gravity action ( |4.1|) does 
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allow series of the anisotropic scaling solutions. To give the logarithmic behavior of the 
holographic entanglement entropy, eq.( |3.18|) needs to be satisfied necessarily. Hence we 
can generally find three solutions for any fixed d (d > 1). Unfortunately, by checking 
it carefully, we haven't found any solution meeting the non-trivial conditions ( |3.20| )- 
( |3.24|) . Thus, it seems that the gravity background ( |3.1| ) which has anisotropic scaling 
property in deep IR with the logarithmic behavior of the holographic entanglement en- 
tropy haven't been covered by above constructions. One possible remedial measure is 
when the anisotropy is only perturbations on the isotropic background, the IR solutions 
(|3.15|) will no longer valid. Instead a(r) behaves like eq. (|3.26|) in the full geometry. Since 
a"{r)/a{r) ~ a'{r)/a{r) ~ 5{a{r)) / a{r) <C 1, we find 0(r) ~ Vt{r) ~ 5{a{r)) / a{r) <C 1. 
If ignoring this small quantity in eq. (PID - (P^) , especially Z^ 1 = which leads H 
field decoupled, the action ( |4. 1| ) will reduce to the standard Einstein-Maxwell-dilaton 
theory and so does the solutions of EOM! One can actually verify that the dilaton 
equation is now satisfied automatically (as emphasized in Ref.[Q). On the other hand, 
even if keeping the first order quantity, it is also approximately satisfied, with only first 
order corrections at most, not broken seriously. Hence the dilaton constraint is sub- 
stantially reduction. Our effective gravity model does realize systems with perturbative 
anisotropy, including those with Fermi surfaces absolutely. 



5. Conclusion 

In this paper, we study the hidden Fermi surfaces in holography by searching the 
logarithmic behavior of the entanglement entropy, without probe fermions in the bulk, 
which was first proposed in Ref.0. We successfully construct the purely classical 
gravity duals for Fermi and non-Fermi liquids for isotropic and anisotropic systems. In 
both systems, the leading contribution of the finite part of the holographic entanglement 
entropy comes from either deep IR region or the deep UV region of the minimal area 
surface. It has been shown that only part of non- Fermi liquids is allowed for isotropic 
systems when the IR piece is dominant in Ref . || while we explicitly prove that it is not 
true for generically anisotropic systems. When the deep UV region contributes to the 
logarithmic behavior of the entanglement entropy, both Fermi and non-Fermi liquids 
can be constructed in both systems. Furthermore, the hyperscaling violation exponent 
is found to be 9 = d(d — l)/{d — 2), d > 2 in isotropic scaling geometries with Fermi 



surfaces in the UV case instead of 9 — d — 1 in the IR case0, || [19| . In the end, we 
also construct a gravity model for anisotropic background which works effectively for 
perturbative anisotropy and does allow series of anisotropic scaling solutions in deep 
bulk interior. However, the dilaton equation is not satisfied automatically and leads to 
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a constraint on the metric. For the purpose of searching a full anisotropic solution, one 
needs a better gravity model. We leave it to future study. 



6. Appendix 

Let's start from the approximated expression for the entanglement entropy 

2 d - 1 R d L d - 1 f 1 du R — — 

S A « / -Wl + CAu) 2 " (6-1) 

where e is the UV cut-off (e — > in the UV limit), leading to the divergent part of 
the entanglement entropy which agrees with area law. Recall that we simply drop the 
u 2d term in the denominator of eq. fl2.22|) to obtain this formula, with the assumption 



that the UV contribution dominant. Since we are only interested in the leading finite 
part of the entanglement entropy, we will drop the UV cut-off in the lower bound of 
the integral. Instead we introduce a new UV scale parameter defined by ruv = rp/A s , 
where 5 is a positive constant to ensure ruv *C rp. 
Define a new variable 



9 = v 7 ! + (A^) 2n (6.2) 
Substitute eq. (|5.2|) into eq. ( |6.1|) , we deduce 



{S A )fin = —i / d6 77; TTTT^T (6-3) 

nr F J 01 (Q 2 - l) 1+ ^r 



where the integrated bounds are 0! = y/l + ( ruv /r F ) 2n , 6 2 = y/1 + A 2n . When 
(d — l)/2n = 1, we can more explicitly write down the integrated expression as follows 



nd-l TjdTd-l /"fell 1 1 

Now we can readily obtain the result of the integral 

2 d - 1 R d L d - 1 r F 
{Sa) ^= 2r d F ^ g ^ + - (6 " 5) 

where the dots denotes the small quantity of order 0(A°). Evidently, by properly 
choosing the UV scale ruv = rp/A s <C rp, we can exactly derive the logarithmic 
behavior of the entanglement entropy, as it was shown in eq.( |2.24 ). 
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